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Abstract. We study the influence of lepton asymmetry on the evolution of the early 
Universe. The lepton asymmetry I is poorly constrained by observations and might 
^ — , ■ be orders of magnitude larger than the baryon asymmetry b, \l\/b < 2 x 10 s . We 

find that lepton asymmetries that are large compared to the tiny baryon asymmetry 
can influence the dynamics of the QCD phase transition significantly. The cosmic 
i c*j trajectory in the (1b — T phase diagram of strongly interacting matter becomes a 

Oh | function of lepton (flavour) asymmetry. Large lepton asymmetry could lead to a cosmic 

O |' QCD phase transition of first order. 

^ ■ PACS numbers: 98.80Es, 12.38Aw 
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1. Introduction 

O ' 
0\ 

One of the most poorly constrained parameters of cosmology is the specific lepton 
asymmetry I, the lepton number per entropy, mainly because the cosmic neutrino 
background escapes all direct attempts to be measured. The strongest constrains stem 
from a combined analysis of big bang nucleosynthesis (BBN) and the cosmic microwave 
background (CMB), leading to \l\ < 0.02 [HI21E]- 

Cosmological inflation predicts for the observable universe that both, baryon 
and lepton numbers are approximately zero. Therefore baryo- or leptogenesis or a 
combination of both is needed to result in the today measured specific baryon asymmetry 
b ~ 10~ 10 . Before the electroweak transition at Tew ~ 200 GeV, sphaleron processes 
allow the violation of baryon and lepton number in a way, that A(B — L) =0. Within 
the standard model of particle physics (SM) this would lead to I = —(51/28)6, where the 
numerical coefficient depends on the particle content of the universe before the onset 
of the electroweak transition [4J. But it may also be that the asymmetries differ by 
orders of magnitudes. Late (compared to the epoch when sphaleron processes cease to 
be efficient) baryogenesis would lead to b ^> |/|, while late leptogenesis would lead to 
|/| 3> b [HI El [71 El [9] . From the experimental point of view, such a large specific lepton 
asymmetry may be hidden in todays neutrino sector. 
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An likewise unknown, although most dramatic event in the early universe is 
the quark-hadron transition at approximately 10/is after the big bang, where quarks 
and gluons condense to hadrons. The SM predicts a spontaneous breaking of the 
chiral symmetry of QCD and a confining of quarks into hadrons at a (pseudo-) critical 
Temperature T c . This epoch is one of the most interesting in the early universe, but 
experimentally only little is known about this QCD transition at T c w 190 MeV [HI [15]. 
It is well established that the order of the transition depends on the baryon density, or 
equivalently the baryochemical potential which is the energy needed to add one 
unit of baryon number to the system. It is common to draw a QCD phase diagram in 
the ijlb — T plane, in which the cosmic QCD transition is commonly assumed to take 
place at T = T c and /xg ~ 0. We show in this work that this assumption holds only if 
|/| = 0(6). 

The QCD transition is in the focus of the relativistic heavy- ion research programs 
at RHIC and LHC. These experiments are taking place at very small baryochemical 
potential and high temperatures, presumably far away from a critical point. More 
information about a possible first-order transition at large values of \i& are expected 
from the future FAIR program. Todays ordinary matter is located at T « and [is 
equal to the nucleon mass. Extremely dense matter, like the interior of compact stars 
with T < 10 MeV, are expected to reach some colour superconducting phase [TO] . 

For T and \xb being of the order of the QCD confining scale Aqcd, the most reliable 
tool is lattice simulations with quarks. But as a tribute to different possible solutions 
to the so called 'sign problem' [121 02], different phase diagrams with different T c and 
different critical lines have been proposed [UJ [T2J [18] . The widely expected shape of the 
QCD phase diagram in a \Lb — T plane as a combination of results from nuclear physics, 
perturbative calculations and lattice simulations is shown in figure [TJ 

For cosmology, gaining better knowledge of the QCD transition would be great. 
The transition sets the initial condition for big bang nucleosynthesis (BBN). Standard 
BBN relies on the assumption that the Universe can be treated as a homogeneous and 
isentropic radiation fluid at the moment of BBN. Therefore it is of great interest to know 
exactly what happens at the QCD transition [21J. The cosmic QCD transition might 
also create some relics observable today. Most of them can only be generated during 
a first-order transition, like quark nuggets or magnetic fields [23[ [29] , while QCD balls 
could be generated independently of the order of the transition [23]. The transition 
might also lead to the formation of cold dark matter (CDM) clumps [25] and definitely 
leads to a modification of the primordial background of gravitational waves [26], [28] . 

In this work, we trace the baryo-, lepto- and charge chemical potentials in the early 
Universe as a function of temperature around the QCD phase transition. We assume 
that the Universe is neutral and that electric charge, baryon and lepton (flavor) numbers 
are conserved. For vanishing lepton asymmetry this situation has also been studied in 
[33] and a similar work can be found in |32j . 

In section 2 we explain our assumptions about the particle content and their 
interactions in the early universe after the electroweak transition and before BBN. 
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Figure 1. Sketch of the phase diagram for physical QCD, based on findings from 
nuclear physics, lattice QCD and perturbative calculations. Solid lines indicate a first- 
order phase transition, while the dashed line indicates a crossover. The exact phase 
diagram with or without a critical end point (T e ,fj, e ) is still under debate. So is also 
the (pseudo-)critical temperature T c = 192 (7) (4) MeV [IH [15], which differs from 
T c = 164 ± 2 MeV, the value found in [16]. While [18] argue that it is still an open 
question if a critical end point exists at all, [19] find /i e = 360±40 MeV and T c — 162±2 
MeV. The calculations and methods leading to these different results are discussed and 
compared in [HI [H [H [12] . 

We introduce the relevant conserved quantum numbers and use the corresponding 
conservation laws to solve for the net particle densities and corresponding particle 
chemical potentials in the quark and hadron phases. In the third section, we introduce 
the chemical potentials related to the conserved quantum numbers and study the 
trajectories of the Universe in the corresponding \i—T planes, depending on the assumed 
lepton asymmetry. We discus some consequences of our findings on the cosmic QCD 
transition in section 4 and conclude in section 5. 

2. Evolution of particle chemical potentials and densities 

In the early Universe, after the electroweak phase transition at Tew ~ 200 GeV and 
before the onset of neutrino oscillations at a few MeV, baryon number B, lepton flavour 
numbers Lf, and the electric charge Q are conserved. Additionally, we assume that 
Q = 0. A neutral universe seems to be reasonable, because several observations point 
to neutrality and lack of currents on large scales. The work of [30] shows that a possible 
charge asymmetry is annihilated in less than a Hubble time for 100 GeV > T > 1 eV. 
It is shown too, that currents are damped for all temperatures T > 1 eV. 
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We also assume that all globally conserved quantum numbers are also conserved 
locally. This means that our approximation is applicable at length scales larger than the 
largest scale on which transport phenomena can show up, given we apply homogeneous 
initial conditions. The largest mean free path is that of neutrinos and thus the local 
physics on scales less than the neutrino mean free path might differ from the results 
obtained in this work. 

The net particle density n,i is the density of a particle minus the density of its anti- 
particle. In thermal and chemical equilibrium (and neglecting effects of interactions) it 
can be described as 

m(T,u i ) = — [ EJE 2 - mi 2 ( =^ =rj | dE, (1) 

1 ^ ^J mi V l v exp(^)±l exp(^)±i; KJ 

for a particle with mass m^, chemical potential \X{ and helicity degrees of freedom g^. The 
sign is fixed by the spin of the particle, + for fermions, — for bosons. As the Universe 
expands and cools, the particles annihilate at temperatures comparable to their mass 
and eventually fall out of equilibrium. 

Lifting the global conservation laws to local ones we obtain equations for the specific 
lepton flavour asymmetry, specific baryon asymmetry and the charge density: 

nf + n Uf 

l f = for / = e, //, r, (2a) 

b = with hi = baryon number of species i, (26) 

^-^ s 

i 

= qiUi with qi = charge of species i, (2 c) 

i 

with the entropy density s = s(T). For any given temperature T, the free parameters in 
this set of equations are the chemical potentials in the net particle densities rii(T, //j) and 
the specific lepton flavour asymmetries If. The baryon number b w (8.85 ±0.24) x 1CT 11 
is given, based on [34J. 



2.1. Quark phase 

In the following we describe the evolution of the net particle densities of all SM particles 
in the quark phase (T > Tqcd)- 

In the early Universe the Hubble expansion rate H is always bigger than the particle 
interaction rate V and the particle densities are usually not high enough for many body 
reactions. The main processes that keep all particles in thermal and chemical equilibrium 
are 2 <-> 2-processes. So, all particles in the quark phase and in the hadron phase are 
in local thermal (LTE) and local chemical equilibrium (LCE). The chemical potentials 
of photons and gluons are zero /x 7 = \i g = 0, as they can be produced and annihilated 
in any number. As a consequence of that, the chemical potentials of a particle i and its 
anti-particle i are equal, but of opposite sign: ^ = —ft-. 
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The leptons are then linked via reactions like fa + fj ^ + Uf. for all combinations 
of flavours, leading to relations like 

f*e - fa = - Pup- (3) 

Due to flavour mixing in the quark sector, we only need to distinguish between up 
and down type quarks with chemical potentials \i u = \i c = fi t for the up, charm and 
top quarks and fid = fi s = fib for the down, strange and bottom quarks. These quark 
chemical potentials are further linked via weak reactions like u + e~ f± d + i/ e , which 
leads to fi d = fi u + fi e — fi Ue . 

The number of independent chemical potentials in the quark phase is thus reduced 
to four leptonic potentials (fi e , fi Ue , fi Ufl , fi Vr ) and just one for the quarks (fi u )- For 
temperatures well below raw and in chemical and thermal equilibrium, the evolution of 
the net particle densities as a function of temperature is described by five conservation 
equations 

= n e (T,fx e ) +n Ve (T,fi Ue ) - l e s(T), (4a) 
= n M (T, fi v ^fi e , fi Ue ) + n Vl ,(T,fi vli ) - Z M s(T), (46) 
= nr(T, /JL e , /JL„ e , fj,„ T ) +n UT (T,fi UT ) - l T s(T), (4c) 
= n u (T, fi u ) + n d (T, ji u , ji e , fi Ve ) + n c (T, fi u ) + n s (T, fi u , fi e , fi Ve ) (Ad) 

+ n b (T, fi u , fi e , fi Ue ) - 3bs(T), 
= n u (T, fi u ) + n c (T, fi u ) - n e (T, fi e ) - n M (T, fJ,^, fM e , fJ>v e ) (4e) 

- n r (T, fM e , fj,^, fj,„ T ) - bs(T). 

The last equation combines the conservation of electric charge and baryon number in a 
convenient manner. Exact solutions to this system of coupled integral equations must 
be obtained numerically. However, before doing so, we estimate the evolution of the 
chemical potentials at high temperatures via an analytic approach. 

At T > T c , but close to the QCD transition, we consider the three light quarks u, 
d and s, as well as electrons and myons with their corresponding neutrinos (the masses 
of these particles are neglected). Tau leptons annihilate well before the QCD transition 
(at T T ~ 600 MeV), hence tau neutrinos are decoupled at the temperatures of interest. 
For relativistic particles (T ^> mj,/ij) the net number densities become 

n i — -9iT 2 fti + O(fi^) for bosons, (5) 
3 

= -giT 2 fii + 0(f£) for fermions, (6) 
6 

while we neglect the nonrelativistic particles. With these approximations equations ( |4al) 
to (I4~cj) become: 

l e s(T) = ^(2fi e + fi Ue )T 2 } (7a) 

l,s(T) = ^(2fi, + fi^)T 2 , (7b) 

l T s(T) = \^T 2 . (7c) 
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Let us further assume (for simplicity) that l e = Z M = l T = 1/3, which implies, together 
with ([3]), that electrons and muons, as well as electron and muon neutrinos have equal 
chemical potentials, while (TTcjl fixes fj,„ T = ls(T)/T 2 . Consequently, from the complete 
set of conservation equations, the chemical potentials for all particles can be expressed 
in terms of the lepton and baryon asymmetries, I and b (for three quark and two lepton 
flavours) : 

to lbs(T) lls(T) 



T 2 T 3 6 T 3 

to = 1 bs(T) 1 ls(T) 

T 2 T 3 12 T 3 

fl e 1 yL V e 1 ls { T ) 



(8b) 

T IT 4 T 3 ^ 
While the chemical potentials for the leptons ( |8cj) depend on the lepton asymmetry 
only, the quark chemical potentials given by (I8aj) and (j86|) depend on baryon and lepton 
asymmetry. 

We would like to stress that the result flHoj) to fl8cj) strongly depends on the number 
of relativistic quark and lepton flavours. E.g. increasing the temperature to T ~ 1 
GeV, the inclusion of tau leptons and charm quarks modifies the coefficients in the 
above result, i.e. 

to 9 bs(T) 1 ls(T) 



T 32 T 3 8 T 3 
to 15 bs(T) lls(T) 

~f ~ 32^ 8 T 3 
to Hs{T) 1 bs(T) 
f~ ~ 4 T 3 + 16 T 3 
/i, e 1 Zs(T) 1 6s(T) 



(9 a) 
(96) 
(9c) 
(9d) 



T 2 T 3 8 T 3 

It actually turns out that the previously found independence of the chemical potentials 
of the leptons from the baryon asymmetry is a coincidence of the three quark flavour 
case. The deeper reason is that the charges of one down and one strange quark just 
compensate the charge of one up quark. In general, all particle chemical potentials 
depend on b and I. For three quark flavours this is also the case when the strange quark 
mass is taken into account. 

For physical particle masses (e.g. the strange quark mass and the mass of the muon 
cannot be neglected during the QCD transition) an analytic solution is not possible. 
The full numerical solution is presented below. 



2.2. Hadron phase 

The hadron phase (for T < T c ) contains besides hadrons, electrons and muons as well 
as all three types of neutrinos. At low temperatures all quarks are bound in neutrons 
and protons, but close to the QCD transition also mesons, mainly pions, are produced 
in significant numbers. The conservation of quantum numbers, expressed in terms of 
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the net particle densities depending on four independent chemical potentials, gives rise 
to the relations: 

= n e (T,fi e ) - n Ve (T, (jL Ue ) - l e s{T), (10a) 
= n^(T,fi p ,n n ,n^) - n^(T, //„J - / M s(T), (106) 
= n p (T,/i p ) +n„(T,/i n ) - bs(T), (10c) 
= n p (T, fip) - n e (T, /i e ) - n M (T, /i p , fi n , fi u J - n n (T, fi p , fi n ). (10 d) 

For a simplified analytic treatment, we neglect all mesons and muons. Only protons, 
neutrons, electrons and electron neutrinos are taken into account. The leptons are still 
relativistic (as long as T > m e /3), while protons and neutrons are non-relativistic. Their 
net particle density is approximated as 

/ miT\ 3 ^ 2 ( / -mA 
Ui ~ 2gi \ l^T J \T ) 6XP V ~T r ) ' 

The system of equations fllOap . (110 cp . and (1104 then becomes (as above, we assume 

le = 1/3) 

ls(T) =2T 2 (2// e + /0, (12a) 

fo(T) = 4 (^) 3/2eXP (^) ( Sinh T + Sinh f ) ' (126) 
^(^^exp^smh^jTV, (12c) 

where we assumed m n « m p = m. The weak reactions e + p ^ n + u e , which are in 
equilibrium at energies above ~ 1 MeV, f 1 1 2 a[) implies 

H e _ ls(T) fi n - fip 

T 6T 3 3T 1 ' 

For /ip/T <C 1 and \x n /T <C 1, we may simplify ( 11261) to obtain 



/i p /i n 6s(T) 

— H = — ; — -exp m/i , 14 

T T c(T) l[ 1 J ' v ; 

with c(T) = 4(^) 3 / 2 . Together with (TT3"]) and (112 cp the proton chemical potential reads 
/jp _ ^ + Ifyexp [m/T] 

T 1 + ^) eX p[-m/T]- 1 ] 

Note that this expression holds as long as fi p /T <C 1. Nevertheless, the expression allows 
us to understand why the lepton asymmetry does not couple to the proton and neutron 
chemical potentials. As T decreases, the influence of the baryon asymmetry eventually 
overcomes any lepton asymmetry, as the term cx 6exp[m/T] grows exponentially. Thus, 
in the low temperature regime, large lepton asymmetries play no role for the proton 
chemical potential. In the same limit, fi p ~ fi n . From equation (TLB"]) we find that a large 
lepton asymmetry results in a large neutrino and electron asymmetries, as fi Ue ~ fi e . 
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As the temperature drops further, /x p /T becomes large and the above approximation 
breaks down. In that case we can approximate sinh(/i/T) in (I12&I) by exp(/x/T)/2 and 
we finally find that 

H P ^m-T\a[c{T)/2bs{T)\. (16) 

In the low temperature regime \i p is linear in T and for small T it runs against the mass 
m. The electron asymmetry finally is just the same as the proton asymmetry and a 
possible large lepton asymmetry is turning into a large neutrino asymmetry in the late 
Universe. 

Thus, large lepton asymmetries in the low temperature regime of the hadron gas 
affect only the chemical potentials of leptons and play no role for \i v and fi n . The full 
calculation for all particles of the standard model with physical masses can again only 
be done numerical, as is shown in the following chapter. 



2.3. Numerical results 

Taking all particle masses into account, equations (14 ap to (J4ej) and fllOal) to (llOcjj) have 
to been solved numerically. We use Gauss-Laguerre integrations for the net particle 
densities and a multidimensional secant method, called Broydn's method [36J. bf To 
avoid numerical instabilities the net particle densities in ([1]) are rewritten as 

(Tm) 3/2 I(T,/i), (17) 

where 

(1 + -X)x/(1 + -X)X 

: m ; v l vlI dx (is) 

(exp[x - f ] + exp[-f ]) (exp[x + £] + exp[-f ]) ' V ; 

For the quark phase we consider the temperature interval 10 GeV > T > 10 MeV 
with all SM particles, even if some of them, like the top-quark or the W ± and Z-bosons 
are to heavy to play a significant role in this temperature regime. All particle masses 
are adopted from the Particle Data Group [35]. We restrict our study to the case of 
an equal flavour asymmetries l e = l^ = l T = 1/3, the baryon asymmetry is fixed to 
b = 9 x 10~ n in all our numerical calculations. As already discussed above, if sphaleron 
processes are efficient and no additional baryon or leptogenesis happens after they stop, 
the generic value of I = — [4J. This leads to a universe dominated by antimatter in 
the leptonic sector and matter in the baryonic sector. Below we discuss two examples 
/ = -b and I = 3 x 10" 4 . 

The evolution of the net lepton densities in the quark-gluon phase is shown in 
figure [2] for I = —b. At high temperatures there is no difference between the three 
lepton flavours. A negative net particle density means, that more antiparticles than 
particles exists. Below ~ m T /3 ~ 600 MeV the tau leptons disappear, giving rise to an 
increase of net tau neutrino density, in order to keep l T constant. At the same time the 
number of positrons and positive muons increases in order to balance the positive charge 
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Figure 2. Evolution of net lepton densities n, in the quark phase for I — —b. We plot 
them with respect to the entropy density s(T). 
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Figure 3. Evolution of net quark densities in the quark phase for / = —6. 
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that can no longer reside in the tau sector. To conserve flavour this is compensated by 
a corresponding change in the neutrino sector. 

The net quark densities for this case are shown in figure El The two different quark 
types u and d start at different points. If the temperature reaches the mass threshold 
of a quark species, it annihilates or decays and the remaining up- or down-like quarks 
increase their number densities to keep charge neutrality. For example the strange and 
down quarks react on the bottom disappearance. The picture shows, that considering 
only u, d, and s quarks around the QCD transition at Tqcd ~ 200 MeV is a reasonable 
approximation. 

Increasing the lepton asymmetry leads to drastically changes of net particle 
densities. In figures H] and [5] their evolution in the quark-gluon phase are shown for 
I — 3 X 10 -4 . A positive lepton asymmetry leads to a dominance of negatively charged 
leptons. It comes without surprise that the net lepton densities react in a direct way, 
as they obviously depend on If. However, the quarks change their net densities by an 
comparable amount, as already indicated by our previous estimates (l8ajl and ( l86j) . The 
net down-like quark densities become negative to react on the higher densities of the up- 
type quarks. The densities of the up types has to increase as an reaction of large lepton 
asymmetry to reach the charge neutrality. The different signs of the quark chemical 
potential for large I are also indicated by (I8aj) and (|8foj) . Our numerical calculations 
confirm the analytic approach (|8aj) to (|8cj) and show that large lepton asymmetries 
influence the net quark densities strongly. 

Let us now turn to the hadron phase. The numerical solution for all relevant 
particles in the hadron phase are shown in figures ([6]) and ([7]). We consider the evolution 
form 300 MeV down to 1 MeV. Below T ~ 10 MeV effects of neutrino oscillations become 
relevant, which are not included in this work. For these calculations we took all hadrons 
listed by the Particle Data Group up to the mass of = 1414 MeV into account. 

The figures show the evolution of net particle densities for I = —6 and / = 3 x 10~ 4 , 
respectively. Again it can be seen, that I = —b leads to an antimatter dominated universe 
in the leptonic sector until the pion disappearance (not shown in the figure). At the 
lowest temperatures considered, the electron density is strongly linked with the baryon 
number via the charge neutrality of the Universe. The inclusion of light mesons like 
pions and kaons is also important to understand the split between proton and neutron 
densities at temperatures close to the QCD transition. This is another effect that is 
not obvious from the very beginning. As pions are effectively massless, they can carry 
negative charge density in order to compensate for the positive charge density in protons 
and positrons. The influence of larger I on the net particle densities of the neutron and 
proton can be clearly seen, as the evolution differs in the high temperature regime. For 
low temperatures their evolution is independent of Z, as indicated by (TT5]) . The effect of 
large I on fi p and fi n for small T is negligibly. 

Our numerical calculation fits the analytic approaches for both phases well. There 
is clearly an influence of large I on quarks and hadrons at the QCD transition. In 
the following section we take a closer look at the chemical potentials around the QCD 
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Figure 4. Evolution of net lepton densities in the quark phase for / = 3 x 10~ 4 . The 
order of magnitudes is very different from the case / = —6 and the sign is reversed. 
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Figure 5. Evolution of net quark densities in the quark phase for / = 3 x 10~ 4 . In 
comparison to the I = —b scenario, we find a very different behaviour of the net quark 
densities. Charge neutrality and the tiny baryon asymmetry force the net densities of 
down- like quarks to negative values. 
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Figure 6. Evolution of net particle densities in the hadron phase for I = —b. 
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Figure 7. Evolution of net particle densities in the hadron phase for I — 3 x 10~ 4 . 
The influence of a larger I on the net particle densities of the neutron and proton can 
be seen clearly. 
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transition to see this influence more clearly. 



3. Chemical potentials 

Each conserved quantum number can be associated with a chemical potential. The 
particles contribution to the free energy can then be described as 

fi Q n Q + fj, B n B + ^2 ^L f n Lf T> ^ CD ^ N n q + ^2 M ni + ^ n g 

f 1 l 9 

b m I 

with Q denoting charge, B baryon number, and Lf lepton flavour number on the 
l.h.s. and q quarks, 1 leptons, g massive gauge bosons (the gluon and photon chemical 
potentials vanish), b baryons and m mesons. 

The particle fluid in the quark phase contains all types of quarks and leptons. 
The massive gauge bosons are way too heavy to play a crucial role at the QCD epoch, 
nevertheless we keep them in our model. So we get: 

2 12 12 1 

/Mg n « - + - -n s + -n t - -n h - n e - - n T + n w ) 



+ Hb l n i^j + Yl ^ L f ( n f + Uv f 



) 



= ^2 N n Q + Y ^ lUl + l^wnw 

q I 

The Z-boson does not show up, as at high energies its chemical potential vanishes (like 

the one of photons and gluons) and there are no Z-bosons at low temperatures. For the 

W we pick the convention that n w counts W + with a plus sign. A comparison of the 

coefficients in front of the net particle densities gives: 

2 1 . , 

Hu = He = th = g/^Q + gA*B, (19a) 

H d = n a = fi b = + (196) 

Hf = ~HQ+HL f , (19c) 
Vv f = f*L f , (19fi?) 

Hw = HQ- (19e) 

Consequently, the chemical potentials of the globally conserved quantities in the quark 
phase are 

H B {T > T QCD ) = im + 2fi d , (20a) 
H Q (T > T QCD ) = Hu- Hd, (206) 
HL t (T>T qpD )=n vr (20c) 
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b m I 



As above, a comparison of the coefficients leads to 



li B {T < T QCD ) = /i n , 

^QyT < Tqcd) = Hp — = Htt, 

fj, Lf (T < Tqcd) = Hv r 



(21a) 
(216) 
(21c) 



In the following we take a closer look on the chemical potentials of the conserved 
charges and on the effect of a large lepton asymmetry on them at temperatures close to 
the QCD transition. 

3.1. Baryochemical potential 

The baryochemical potential /ie is the energy needed to add a baryon to a thermalised 
state at fixed volume and entropy. In general, it depends linearly on b and I. 

In the quark phase, taking only three quarks (u,d,s), electrons and muons into 
account (and neglecting their masses), we insert ( |8aj) and ( |8fej) into (120 ah to find 
Hb{T > Tq CD ) = 3bs(T)/2T 2 . Thus, it seems that the baryochemical potential would 
just depend on the baryon asymmetry, as one could naively expect. However, if we go 
to a slightly higher temperature, such that we have to include charm and tau, we find 
from (Eij) and <^Ej that 



In general, fis — Hb{°) 0- Thus, for \l\ b, the order of magnitude of the baryochemical 
potential is set by the lepton asymmetry, and not as naively expected by the baryon 
asymmetry. 

This can be understood by the following consideration: increasing the lepton flavour 
asymmetries, the corresponding chemical potentials are forced to react, see (TTab to (!7cj) . 
On the other hand, electric charge and baryon number must not be changed. A large 
net charge in the leptons must be compensated by a large and opposite net charge 
carried by quarks, in a way that the baryon asymmetry remains tiny. That makes it 
very difficult to add a baryon without disturbing that fine balance, and that is why the 
baryochemical potential (120 qp grows as large as the lepton asymmetry. 

Neglecting the physical masses of the particles would have a strong effect on /xg, 
as shown in figured While the simple analytic approach with m = gives \xb oc T (see 
( f22|) ). the numerical results with physical quark masses shows in the low temperature 
regime of the quark phase a different behaviour. The non-relativistic particles give rise 
to an increase of hb below T ~ 3 GeV, i.e. charm, bottom quarks and rs play an 
important role. However neglecting the heavy quarks, but taking the physical masses of 




(22) 
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Figure 8. Trajectory of the quark phase in the — T plane for the I = —b 
scenario. We compare several approximations to the exact result (black line). At 
high temperatures all particles can be assumed to be massless (blue line) , while below 
~ 5 GeV mass thresholds are important. Below 1 GeV the universe is well described 
by three quark (up, down, strange) and two lepton (electron, muon) flavours only 
(magenta line). 



the strange quark and the muon into account gives a reasonable approximation at the 
QCD scale (see figure [8]). 

In the hadron phase at low temperatures, but above temperatures of a few MeV, the 
difference in the masses and chemical potentials of the proton and neutron is negligible 
and one can assume m = m p « m n and /^b = yU„ ~ as given by (fl6l) after pions and 
muons disappeared from the thermal bath. Thus, the baryochemical potential becomes 

c(T) 



Hb(T < rn n /3) ~ m — Tin 



(23) 



2bs(T) 

In the low temperature regime, \xb is linear in T and approaches the nucleon mass m. 

Let us now turn to a more detailed numerical study We evaluate /xs(T) numerically 
in the quark and hadron phases, i.e. expressions (120 ah and (12 1 ah respectively, for various 
values of I. 

If there is no baryo- or leptogenesis after sphaleron processes cease, the standard 
model predicts a negative value for I. The influence of the sign of I on the baryochemical 
potential is displayed in figure [9] In order to illustrate the effect of the sign, we choose 
a rather "large" value for the lepton asymmetrie / = ±3 x 10~ 2 . For negative I, fis is 
positive for all temperatures. In contrast, for positive values of I, ^b can take negative 
values. For the high temperature regime, this is easily seen from ( |22l) . In the hadron 
phase, fiB approaches the nucleon mass at low temperatures, independently of the value 
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Figure 9. Evolution of the baryochemical potential in the quark and hadron phase 
for large lepton asymmetry I = ±3 x 1CP 2 . Lepton asymmetries of different signs 
lead to qualitatively different trajectories. While the trajectories of the quark and 
hadron phase cross for positive /, they miss each other for negative / in the expected 
temperature region at ~ 200 MeV. 




1e-09 1e-08 1e-07 1e-06 1e-05 1e-04 0.001 0.01 0.1 1 

Figure 10. Evolution of the baryochemical potential for negative lepton asymmetries. 
Above the pion threshold, the trajectories of both phases depend strongly on the lepton 
asymmetry. For lower temperatures the lepton asymmetry is negligible. 



Lepton asymmetry and the cosmic QCD transition 



17 



of I. While a positive / always allows for a point of coexistence of the quark and hadron 
phase at a temperature of ~ 200 MeV, the trajectories of the I < case would meet at 
a much lower temperature. 

Figure [10] shows the case of negative lepton asymmetry for different values of I. The 
influence of the lepton flavour asymmetry on the baryochemical potential is evident. The 
point of coexistence moves to larger baryochemical potentials and lower temperatures 
as \l\ increases (not all crossings are shown in the figure). Shortly after the transition 
to the hadron phase, the trajectories are still heavily affected by the lepton asymmetry, 
as pions and other mesons can carry a net charge density. For temperatures below the 
pion threshold, T < m n /3 ~ 50MeV, the hadron phase follows a unique trajectory, as 
predicted by fl23|) . 



3.2. Charge chemical potential 

We associate with the globally conserved and vanishing charge of the Universe the charge 
chemical potential fiQ, which is the energy needed to add a charge unit to a thermalised 
state at fixed volume and entropy. 

For T > Tqcd it depends, like the baryochemical potential, on /i„ and /i^. But 
the effect of a large lepton asymmetry on the evolution is bigger because of the missing 
factor of two for /id- From the approximations (IHaj) and ( |8fej) . it follows for the case of 
three quark and two lepton flavours, fiQ = ls(T)/4T 2 . Taking, as above, the charm and 
tau into account and using ( !9~aj) and ( |96l) we find: 



/i Q (T«r QCD )=(-^ + ^^. (24) 

The charge chemical potential is now depending on b too, and does not even vanish in 
the case b — I, as one could naively expect. Taking the particle masses into account has 
the same effect, namely hq = fiQ^b, I) in general. 

Actually, the lattice simulations that have been used to conclude that the cosmic 
QCD transition would be a crossover have been done with three dynamical quarks. 
They implicitly assume I = and are typically performed at \iq = 0. However, as we 
have seen above, including the mass of the strange quark already spoils that argument. 
We thus have to revisit the question of the order of the cosmic QCD transition. 

After the QCD transition, a charge unit can be added to the Universe in the form 
of either a charged pion, muon or electron. While adding an electron or muon also 
affects the lepton asymmetry, the pions do not care about baryon and lepton number. 
At T « m n /3 ~ 50 MeV the pions are annihilated and the charge neutrality has to be 
ensured by protons and electrons. Using ffT4|) . (fl5l) and (12 1 £>[) one gets 

s(T) 

^ (r<r ^ ) = ( - 96 + %r> + ^e4[- m /T] ^ 
The dependence on the sign of I is shown in figure [TTJ In the quark phase it behaves 
like hq oc l/T 2 , so that the sign of I fixes the sign of [1q. Figure [121 shows the dependence 
of hq for negative I. A negative value of / implies a negative charge chemical potential. 
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Figure 11. Evolution of the charge chemical potential for I = ±3 x 10~ 2 . The 
trajectories of the quark and hadron phase miss each other for both signs of I. 
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Figure 12. Evolution of the charge chemical potential for negative values of I. 



It is interesting that there does not exist any point of coexistence for all the studied 
cases. 
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Figure 13. Evolution of the leptochemical potential associated with the first lepton 
generation for I = ±3 x 10~ 2 . 



3. 3. Leptochemical potentials 

With each lepton flavour a chemical potential is associated. It is the energy needed to 
add a lepton of a certain flavour to a thermalised state at fixed volume and entropy. 
The cheapest way to add a lepton of some flavour is to add a corresponding neutrino, 
see (120 d) and ( 12 left . While for the case of three light quarks and two charged leptons 
at temperatures just above the QCD transition, [i Le = \i L = /Ul t /2 = ls(T)/2T 2 from 
(|8d) and (ITcj) . including the charm and tau at higher temperatures leads to ( |9dl) 

MT«T QCD )= {l~\h) S -^ (26) 

for all three lepton flavours. Again, taking the finite masses of quarks and charged 
lepton into account shows that fiL f = HL f {b,l). 

After the annihilation of the tau and myon, the flavour asymmetries survives in 
the neutrino sector. Only the asymmetry in the electron flavour can be divided into 
electron and electron neutrino asymmetry. This changes with the onset of the neutrino 
oscillations at T ~ lOMeV. The flavours oscillate and the three neutrino sectors become 
linked via v e ^ f± v t and lepton flavour is no longer a conserved quantum number. 
However, the total lepton number still is (if neutrinos are not their own anti-particles, 
i.e. are not Majorana particles). 
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Figure 14. Evolution of the leptochemical potential associated with the first lepton 
generation for negative I. 

4. Consequences for the cosmic QCD transition 

Our findings question some of the established results about the dynamics of the cosmic 
QCD transition. For the scenario I ~ — b we found for the baryochemical potential 
a crossing of the quark and hadron phase trajectories at T m 190 MeV (figure [TO]) . 
This seems to be consistent with a crossover found by lattice simulations, but the 
trajectories for the charge and leptochemical potentials do not cross at T ~ 200 MeV. 
For a homogeneous universe this would require a jump from the quark- to the hadron 
trajectory, which would lead to a phase transition. In our simple model, due to /xq 7^ 
a phase transition for / ~ —b cannot be excluded. However, as we have neglected 
interactions in our analysis, the trajectories might be modified close to T ~ Aqcd and 
our findings just tell us that the situation of /jq 7^ and jii f 7^ has to be studied 
carefully. For the tiny baryon asymmetry and a small lepton asymmetry it would in 
principle be possible to revisit the question of the order of the QCD transition at \iq 7^ 
by means of lattice QCD. 

So far it seems, that increasing the lepton asymmetry leads to a decrease of the 
critical temperature. The opposite effect occurs at the electroweak transition, shown in 
[37]. 

The sign of a possible I asymmetry is important for the dynamics of the transition, 
as is obvious form figures [9l [UJ and [131 However, for large \l\ or for positive I, the 
trajectories of the quark and hadron phases do not cross in the /ig — T plane (at least 
not at T > 100 MeV. This might again suggest that the transition would become first 
order. This is also consistent with the speculation on the existence of a critical end 



Lepton asymmetry and the cosmic QCD transition 



21 



point. Lepton asymmetries of the order of 10 -2 might put the early Universe in the 
vicinity of this critical end point. 

To go into more details is beyond the scope of this work, but it is clear that the 
phase diagram for the cosmic phase transition lives in 5 dimensions for a charge neutral 
universe (3 lepton flavour asymmetries, the baryon asymmetry and temperature). So 
far just a two dimensional slice of it has been studied in greater detail. 

5. Conclusions 

We know very little about the actual lepton asymmetry of the Universe. A large lepton 
asymmetry |/| ^> b is compatible with constraints form primordial nucleosynthesis and 
the CMB if |/| < 0.02. We have shown that such an asymmetry may significantly affect 
the dynamics of the cosmic QCD transition. The cosmic trajectory in the fiB-T plane 
depends on the lepton flavour asymmetries //, besides the baryon asymmetry b. 

Depending on the, yet unknown, structure of the QCD phase diagram and especially 
on the position of a hypothetical critical end point, a large lepton asymmetry might 
result in a first order QCD transition in the early Universe. This seems possible for 
|/| ~ 0.02. Many formerly suggested cosmological consequences of a first order QCD 
transition would be possible (formation of relics, effects on nucleosynthesis, generation 
of magnetic fields, generation of gravitational waves, etc. [29|). 

So far, the conclusions from lattice QCD simulations at fiQ = and /xb <C T have 
been used to conclude that the cosmic QCD transition is a crossover [T7]. We have shown 
that the cosmic trajectory has hq ^ 0. In the case of efficient sphaleron processes hq 
is tiny, in the case of a large lepton asymmetry, fig can be large and it is unclear if the 
conclusion that the cosmic QCD transition is a crossover remains true. For \l\ > 10~ 4 , 
the charge chemical potential differs between the quark and the hadron phase by more 
than 100 MeV at temperatures around the QCD scale (see figure [T2l) . 

A detailed understanding of the consequences of \l\ b on the QCD transition 
could allow us to rule out or find evidence for leptogenesis scenarios that lead to a large 
lepton asymmetry. This is of interest to all scenarios in which the lepton number is 
changed after the sphaleron processes stop to be efficient (late leptogenesis e.g. [9]). 

In this work we assumed that all globally conserved quantum numbers are also 
conserved locally and we have put the focus on the case of equal lepton flavour 
asymmetries. Dropping the latter, obviously leads to a somewhat richer phenomenology, 
which is beyond the scope of this work. While the local conservation of electric charge 
and baryon number are probably excellent approximations, the local conservation of 
lepton number is not correct if distance scales below the neutrino mean free path are 
considered. Consequently, as soon as inhomogeneities become important (and they will 
in the case of a first order phase transitions, as cold spots are more likely to nucleate 
bubbles of the new phase [27]) it would be more realistic to describe an inhomogeneous 
universe with regions of different If. These different regions would be equilibrated via 
the neutrinos, so that possible bubbles of different // at the QCD phase transition 
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would have a minimal radius ^_ mfp ~ 1 m [27]. This has to be compared to the size 
of the Universe at this time, du ~ 10 km [29]. Although giving rise to small effects 
only, precise measurements of the abundance of primordial light elements might be 
sensitive to inhomogeneities produced during the QCD transition, as they might lead 
to inhomogeneous nucleosynthesis. 

Obviously, the consequences of a large lepton asymmetry on the physics of the early 
Universe in the QCD epoch have been overlooked so far. This calls for a reinvestigation 
of many of the suggested effects of the cosmic QCD transition and might allow us 
to improve the limits on the universal lepton flavour asymmetries (before neutrino 
oscillations start). 
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